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Abstract 

Using a sample of 3 x 10^ '0(25') decays recorded by the CLEO detector, we study three-body 
decays of the XcO, Xcij and Xc2 produced in radiative decays of the il){2S). We consider the final 
states 7r"^7r~77, K^K~r}, ppq, tt+tt^t/', K'^K~Tr^, pp'K^, tt'^K'K^, and K'^pA, measuring branching 
fractions or placing upper limits. For Xci 7r+7r~r7, K'^K'tt'^, and 'K'^K~Kg our observed samples 
are large enough to indicate the largest contributions to the substructure. 
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I. INTRODUCTION 



Decays of the Xco, Xa, and Xc2 states are not as well studied experimentally and theoret- 
ically as those of other charmonium states. These states have charge-conjugation eigenvalue 
C = +1, in contrast to the better-studied C = — 1 states J/tp and ip{2S). Their decay 
products thus will differ from those of J ftp and iIj{2S), and may provide complementary 
information on states containing light quarks and gluons. It is possible that the color-octet 
mechanism, ccg 2{qq), could have large effects on the observed decay pattern of the Xcj 
states (H- Assuming that Xcj are the ^Pj cc bound states one would expect that Xco and Xc2 
with J^'^ quantum numbers 0^^ and 2+"*" decay to light quarks via two gluons [2|. Measure- 
ment of any possible Xcj hadronic decays provides valuable information on possible glueball 
dynamics. Thus knowledge of any hadronic decay channels for these states is valuable. 

CLEO has gathered a large sample of e+e~ i^i^S) events, which leads to copious 
production of the Xcj states in radiative decays of the tp{2S). The ip{2S) branching fractions 
have been recently measured with high precision: 

B{^{2S) ^ 7Xco) = (9.22 ± 0.11 ± 0.46)%; (1) 

B{ij{2S) ^ 7Xci) = (9.07 ± 0.11 ± 0.54)%; (2) 

i3(V'(25) ^ 7Xc2) = (9.33 ± 0.14 ± 0.61)%. (3) 

We describe a study of selected three-body hadronic decay modes of the Xcj to two charged 
and one neutral hadron. This is not an exhaustive study of Xcj hadronic decays; we do not 
even comprehensively cover all possible h~^h~h^ decays, but simply take a first look at the 
rich structure of Xcj decays in our initial ipi2S) data sets. A subset of these modes has been 
investigated by BES [4] . With the CLEO III detector configuration [HI , we have recorded an 
integrated luminosity of 2.57 pb~^ and the number of ip{2S) events is 1.56 x 10^. With the 
CLEO-c detector configuration [6| we have recorded 2.89 pb~^, and the number of events 
is 1.52 X 10^. The apparent mis-match of luminosities and event totals is due to different 
beam energy spreads for the two data sets. 

II. MEASUREMENT OF THE BRANCHING FRACTIONS 

Our basic technique is an exclusive whole-event analysis searching for iIj{2S) — »• 7Xcj 
followed by a three body decay of the Xcj to tt^tt't], K^K~ri, ppr], tt+tt"?]', K^K~tt^, 
ppTT^, Ti^K^K^, or K^pA. A photon candidate is combined with three hadrons and their 
4-momentum sum constrained to the known total beam energy and the initial momentum 
caused by the two beams crossing at a small angle taking into account the measured errors 
on the reconstructed charged tracks, neutral hadron, and transition photon. We cut on 
the of this fit, which has four degrees of freedom, as it strongly discriminates between 
background and signal. For most modes we select events with an event 4-momentum fit 
less than 25, but background from i'{2S) J/ipir'^n^ followed by charged two-body decays 
of the J/ip, with one of the 7r° decay photons lost, fakes 4'{2S) — > 7Xcj •yppn^. For this 
mode the cut on is tightened to 12. The efficiency of this cut is ~ 95% for all modes 
except ppTi'^, where it is ~ 80%. 

Efficiencies and backgrounds are studied in a GEANT-based simulation [3| of the detector 
response to e+e~ 4'{2S) events. Our simulated sample is roughly ten times our data 
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sample. The radiated photon is generated according to an angular distribution of 1 + A cos^ 6, 
where 9 is the radiated photon angle relative to the positron beam axis. An El transition, 
as expected for ip{2S) 'JXcj, implies A = 1,-1/3, +1/13 for J = 0,1,2 particles. The 
efficiencies we quote use this simulation, and differ from efficiencies using A = by up to a 
few percent. 

Photon candidates are selected by their energy depositions in the Csl crystal calorimeter. 
They have a transverse shape consistent with that expected for an electromagnetic shower 
without a charged track pointing toward it. They are required to have an energy of at least 
30 MeV. Photon candidates that are used to make neutral particles further must have an 
energy of more than 50 MeV if they are not in the barrel, | cos6^\ > 0.82, of our calorimeter. 
The 7r° — 77 and rj —>■ candidates are formed from two-photon candidates that are 
kinematically fit to the known resonance masses using the event vertex position, determined 
using charged tracks constrained to the beam spot. We select events with a from the 
kinematic mass fit with one degree of freedom of less than 10. Transition photon candidates 
are vetoed if they form a 7r° or 77 candidate when paired with a second photon candidate. 

We also reconstruct the rj 7r+7r~7r° mode combining two charged pions with a vr" — >■ 77 
candidate, increasing the number of 77 candidates by about 25%. The same sort of kinematic 
mass fit as used for vr^'s and 77 — 77 is applied to this mode, and again we select those giving 
a of less than 10. Similarly we combine the mass-constrained rj candidates together with 
two charged pions to make rj' candidates, mass-constrain them, and select those with < 10. 
In addition, we include the decay mode rj' 7p. Here the background is potentially high 
because of the large number of noise photons, so we require -Ephoton > 200 MeV. In addition, 
we require the 7r~^7T~ mass to be within 100 MeV/c^ of the mean p mass. 

Charged tracks satisfy standard requirements [sl that they be of good fit quality. Those 
coming from the origin must have an impact parameter with respect to the beam spot less 
than the greater of (5.0 — 3.8 -p) mm and 1.2 mm, where p is the measured track momentum 
in GeV/c. The tt+tt^ and A —>■ pn^ candidates are formed from good-quality tracks 

that are constrained to come from a common vertex. The Kg flight path is required to be 
greater than 5 mm and the A flight path greater than 3 mm. The mass cut around the 
Kg mass is iblO MeV/c^, and around the A mass ±5 MeV/c^, both about three times the 
resolution. Events with only the exact number of selected tracks are accepted. This selection 
is very efficient, >99.9%, for events passing all other requirements. 

Pions are required to have specific ionization, dE/dx, in the main drift chamber within 
four standard deviations of the expected value for a real pion at the measured momentum. 
For kaons and protons, a combined dE/dx and RICH (ring imaging Cherenkov counter) 
likelihood is formed and kaons are required to be more kaon-like than pion- or proton-like, 
and similarly for protons. Cross feed between hadron species is negligible after all other 
requirements. 

In modes comprising only two charged particles, there are some extra cuts to eliminate 
QED background which produce charged leptons in the final state. Events are rejected if the 
sum over all the charged tracks produces a penetration into the muon system of more than 
five nuclear interaction lengths. Events are rejected if any track has 0.92 < E/p < 1.05 and 
it has a dE/dx consistent with an electron. This latter cut is not used for pp modes because 
anti-protons tend to deposit all their energy in the calorimeter. These cuts are essentially 
100% efficient for the signal, and ensure this QED background is negligible. 

The efficiencies averaged over the CLEO HI and CLEO-c data sets for each mode in- 
cluding the branching fractions rj —>■ 77, 7] 7r"^7r~7r*^, and rj' — * rjTC^TT' [9!] are given in 
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FIG. 1: Mass distribution for candidate 
Xcj 7r+7r~?7 events. The displayed fit is 
described in the text. 
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FIG. 2: Mass distribution for candidate 
Xcj K^K^rj events. The displayed fit is 
described in the text. 
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FIG. 3: Mass distribution for candidate 
Xcj — > Ppf] events. The displayed fit is de- 
scribed in the text. 



FIG. 4: Mass distribution for candidate 
Xcj TT^TT^rj' events. The displayed fit is 
described in the text. 



Tables Ulllnl for Xco, Xa, and Xc2 respectively. Figures [T]l8] show the mass distributions for 
the eight Xcj decay modes selected by the analysis described above. Signals are evident in all 
three Xcj states, but not in all the modes. Backgrounds are small. The mass distributions 
are fit to three signal shapes, Breit-Wigners convolved with Gaussian detector resolutions, 
and a linear background. The Xcj niasses and intrinsic widths are fixed at the values from 
the Particle Data Group compilation The detector resolution is taken from the simula- 
tion discussed above. The simulation properly takes into account the amount of data in the 
two detector configurations, and the distribution of different decay modes we have observed. 
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FIG. 5: Mass distribution for candidate 
Xcj K^K~t:^ events. The displayed fit 
is described in the text. 
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FIG. 6: Mass distribution for candidate 
Xcj ppTT^ events. The displayed fit is de- 
scribed in the text. 
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FIG. 7: Mass distribution for candidate 
Xcj — >■ KgK~Tr~^ events. The displayed fit 
is described in the text. 



FIG. 8: Mass distribution for candidate 
Xcj K~^pA events. The displayed fit is 
described in the text. 



We approximate the resolution with a single Gaussian distribution, and variations are con- 
sidered in the determination of systematic uncertainty. The detector resolution dominates 
for the Xci and Xc2, but is similar to the intrinsic width of the Xco- The fits are displayed in 
Figures [T]l8] and summarized in Tables Ul lllli Note that for the Xco in Tabled the five modes 
for which no significant signal is found (yield less than three standard deviations from zero) 
are forbidden by parity conservation. 

We consider various sources of systematic uncertainties on the yields. We varied the fitting 
procedure by allowing the Xcj niasses and intrinsic widths to float. The fitted masses and 
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TABLE I: Parameters used and results of fits to the Xcj mass distributions of Figures [T]15] for the 
XcO signal. The fit is described in the text and the yield error is only statistical. If no significant 
signal is observed we also show an upper limit at 90% confidence level. 



Mode 


Efficiency (%) 


Resolution (MeV) 


Yield 




18.2 


6.23 


4.2 ±4.1 (< 10.6) 


K+K-rj 


13.8 


6.10 


3.1 ±2.7 (< 9.3) 


pprj 


16.1 


5.42 


1 7 7+5.2 
-^'•'-4.8 




8.0 


4.38 


2.0 ±3.8 (< 8.5) 




27.6 


6.47 


-3.3 ±2.7 (< 4.7) 


ppTT^ 


27.8 


5.80 


46.4t^;^ 




19.8 


4.75 


0.0 ± 1.0 (< 2.7) 


K+pA 


16.8 


4.38 





TABLE IL Parameters used and results of fits to the Xcj mass distributions of Figures [T]l8] for the 
Xci signal. The fit is described in the text and the yield error is only statistical. If no significant 
signal is observed we also show an upper limit at 90% confidence level. 



Mode 


Efficiency (%) 


Resolution (MeV) 


Yield 




18.8 


5.85 


255lig 


K+K-ri 


14.6 


5.66 


14.111^ 


pprj 


17.7 


5.41 


3.2 ±2.3 (< 7.6) 




8.5 


4.37 


57.6l?;4 




29.2 


6.79 


157 ± 13 


ppTT^ 


30.1 


5.23 


9.9;|^ 




20.6 


4.37 


249 ± 16 


K+pA 


17.7 


4.38 


16.3lt^ 



widths agree with the values from the Particle Data Group [9||, and we take the maximum 
variation in the observed yields, ±4%, as a systematic uncertainty from the fit procedure. 
Allowing a curvature term to the background has a negligible effect. For modes with large 
yields we can break up the sample into CLEO III and CLEO-c data sets, and fit with reso- 
lutions and efficiencies appropriate for the individual data sets. We note that the separate 
data sets give consistent efficiency-corrected yields and the summed yield differs by 2% from 
the standard procedure, which is small compared to the ±8% statistical uncertainty. We 
take this as the systematic uncertainty from our resolution model. From studies of other 
processes we assign a ±0.7% uncertainty for the efficiency of finding each charged track, 
±4.0% for the 77 resonances, ±2.0% for each extra photon, ±1.3% for the particle identifi- 
cation for each K and p, ±2.0% for secondary vertex finding, and ±3.0% from the statistical 
uncertainty on the efficiency determined from the simulation. We study the cut on the of 
the event 4-momentum kinematic fit in the three large yield Xci signals by removing the 
cut, selecting events around the Xci mass peak, subtracting a low-mass side band, the only 
one available, and comparing the simulated distribution for signal events with the data 
distribution. This comparison is shown in Figure O The agreement between the data and 
simulation is good, and comparing the inefficiency introduced by our cut on the 4-momentum 
kinematic fit between the data and the simulation we assign a ±3.5% uncertainty on the 
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TABLE III: Parameters used and results of fits to the Xcj mass distributions of Figures [T]15] for the 
Xc2 signal. The fit is described in the text and the yield error is only statistical. If no significant 
signal is observed we also show an upper limit at 90% confidence level. 



Mode 


Efficiency (%) 


Resolution (MeV) 


Yield 


PpT] 
ppTT^ 

TT+K-K2, 
K+pA 


18.5 

14.6 

17.2 

8.5 

27.5 

29.3 

20.2 

17.5 


5.58 
5.53 
5.08 
4.32 
6.85 
5.10 
4.45 
4.32 


26.2+^:^ 

6.9 ± 2.9 (< 12.5) 

^•"3-3.0 

24.8+1? 
37.1^11 
36.8^1^ 
42.1^^1 
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FIG. 9: Distribution of the for the event 4-momentum kinematic fit shown here with a Xci mass 
cut and sideband subtraction. Plot (a) is the 7r"*"7r~r/ mode, (b) is the K'^K^tt^ mode, and (c) is 
the 7r~^K~KQ^ mode. The data are shown by points and the simulation of signal events is shown 
by the solid line normalized to have the same area. 



efficiency due to uncertainty in modeling this distribution. The simulation was generated 
assuming three-body phase-space for the Xcj decay products. Deviations from this are to 
be expected. Based on the results of the Dalitz plot analyses discussed below we correct the 
efficiency in the Xci vr+vr"?], Xci K^K'tt^, and Xci K~n~^Kg modes by a relative 
—2.2%, —1.5%, and +6% respectively to account for the change in the efficiency caused by 
the deviation from a uniform phase space distribution of decay products to what we actually 
observe. We apply an additional ±5% uncertainty on all other modes to account for the 
effect of resonant sub-structure on the efficiency. The total systematic uncertainty is mode 
dependent, but is roughly 10%; it is higher for modes with many photons in the final state 
and lower for those with only one. 

To calculate Xcj branching fractions, we use previous CLEO measurements for the 
ipi'^S) IXcJ branching fractions from Equations [UlS] [3]. The uncertainties on these 
branching fractions are included in the systematic uncertainty on the Xcj branching frac- 
tions we report. Also there is a 3% uncertainty on the number of 'ipi^S) produced. Results 
for the three-body branching fractions are shown in Table IIVI Where the yields do not show 
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TABLE IV: Branching fractions in units of 10^'^. Uncertainties are statistical, systematic due to 
detector effects plus analysis methods, and a separate systematic due to uncertainties in the ip{2S) 
branching fractions. Limits are at the 90% confidence level. 



Mode 


XcO 


Xcl 


Xc2 




< 0.21 


5.0 ±0.3 ±0.4 ±0.3 


0.49 ±0.12 ±0.05 ±0.03 


K+K-r] 


< 0.24 


0.34 ±0.10 ±0.03 ±0.02 


< 0.33 


ppr] 


0.39 ±0.11 ±0.04 ±0.02 


< 0.16 


0.19 ±0.07 ±0.02 ±0.01 




< 0.38 


2.4 ±0.4 ±0.2 ±0.2 


0.51 ±0.18 ±0.05 ±0.03 


K+K--K^ 


< 0.06 


1.95 ±0.16 ±0.18 ±0.14 


0.31 ±0.07 ±0.03 ±0.02 




0.59 ±0.10 ±0.07 ±0.03 


0.12 ±0.05 ±0.01 ±0.01 


0.44 ± 0.08 ± 0.04 ± 0.03 


TT+K-K^ 


< 0.10 


8.1 ±0.6 ±0.6 ±0.5 


1.3 ±0.2 ±0.1 ±0.1 


K+pA 


1.07 ±0.17 ±0.10 ±0.06 


0.33 ±0.09 ±0.03 ±0.02 


0.85 ±0.14 ±0.08 ±0.06 



clear signals we calculate 90% confidence level upper limits using the yield central values 
with the statistical errors from the yield fits combined in quadrature with the systematic 
uncertainties on the efficiencies and other branching fractions. We assume the uncertainty 
is distributed as a Gaussian and the upper limit is the branching fraction value at which 
90% of the integrated area of the Gaussian falls below. We exclude the unphysical region, 
negative branching fractions, for this upper limit calculation. We note that the ratio of rates 
obtained from isospin symmetry (see Appendix IVll Equations [231 and [3T1) . expected to be 
4.0, is consistent with our measurement: 



r(x.i - ir+ir-vrO) -4.2 ±0.7. (4J 

Our results are consistent with branching fractions and upper limits in the it+K~K^ and 
7r+7r~?7 modes from BES [^], but more precise. 



III. SUBSTRUCTURE ANALYSIS 

We perform a Dalitz plot analysis on the modes with the highest statistics, Xci 7r"^7r~?7, 
Xcl — > K'^K'Ti^, and Xci 7i~^K~Kg, to study the two-body substructure. For the Dalitz 
analysis only those events within 10 MeV, roughly two standard deviations, of the observed 
Xcl signal peak mean in the specific mode are accepted. For Xci ^ tt+tt^?] there are 228 
events in this region and the signal fit finds 224.2 signal events and 5.1 combinatorial back- 
ground. For Xcl K~^K~7T^ there are 137 events accepted with the fit finding 137.8 signal 
and 2.4 background events, and for Xci 7i+K~K^, the numbers are 234 events, of which 
233.2 are signal and 0.8 are background. In all cases the contribution from the tail of the 
Xc2 is less than one event. 

An unbinned maximum likelihood fit is used in order to perform the Dalitz plot analysis 
ll| . In order to assess the fit quality we use an adaptive binning technique and calculate 



a probability for Pearson statistics. Efficiencies are determined with simulated event samples 
for the Xcl decay generated uniformly in phase space, and run through the analysis procedure 
described above. The efficiency across the Dalitz plots is fit to a two-dimensional polynomial 
of third order in the Dalitz plot variables. The fits are of good-quality and the efficiency is 
generally flat across the Dalitz plot. 
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When fitting the data small contributions from backgrounds are neglected. We are ex- 
amining the e~^e~ — >• ipi^S) iXci process. In such a decay the Xci should be polarized. In 
principle a complete analysis would take into account the angle of the photon with respect 
to the e"'"e~ beams' collision axis and decompose the Xci decay into its partial waves. We 
use a simple model to analyze our data sample which is adequate for seeing the largest con- 
tributions to the substructure in our small sample. We take the Dalitz plot matrix element 
\M\^ to be a sum of non-interfering resonances, 

\M\' = Y.\AR?-nl. (5) 

R 

The quantity AR{m) represents the amplitude of each resonance contribution with angular 
distributions Vt\ taken from Ref. [13] as shown in Table |Vl 

TABLE V: Angular distributions from Ref. [lol ] used in the present analysis. The notation foUows 
the original publication. We assume the decay of a particle with spin J ( J = 1 for Xci in our case) 
to the resonance R of spin j and a pseudoscalar particle with relative orbital momentum L. The 
angle 6 is the resonance decay angle with respect to the pseudoscalar particle in the resonance's 
rest frame. The term = 7^ — 1 is a relativistic correction factor, where 7ij is the resonance 
Lorentz factor in the Xci rest frame. 



J - 


-.j + L 


Angular distribution, 17^ 


1- 


-^0+1 


uniform 


1- 


-^1+0 


1 + cos^ e 


1- 


-^1+2 


l + (3 + 4z2)cos2 6l 


1- 


-^2+1 


(1 + + 3cos2 e + 9z2(cos2 9 - 1/3)2] 



Narrow resonances are described with a Breit-Wigner amplitude 

AR{m) = — — — . (6) 

m%-m^ - imTR^ total ("^) 

The aji coefficients are fit parameters giving the amplitudes of the resonance with spin j- 
and mass-dependent width 



mTR^ totaiM = Vi3(i? /) ■ mRrJ—] — , (7) 
Y \PrJ m 



where the resonance mass ttir, width TR^nd branching ratio B{R — > /) into the final state 
/ are taken from previous experiments [9|; p and pr are the decay products' momenta in 
the resonance rest frame and its value at m = uir. For the scalar resonance ao(980) we use 
a Flatte parameterization in the style of the Crystal Barrel Collaboration (l3| . 



ao{980) m _ ^2 _ i[g2^p^^(^rn) + g^j^pRKim)] ' 

where tur, g^jn, and g^^ are the resonance mass and coupling constants, m is the invariant 
mass of the resonance final state, and Pab{^) is a phase space factor for the particular final 
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FIG. 10: (a) Dalitz plot, (b)~(d) projections on the three mass squared combinations for Xci 
TT+vr^r/. The displayed fit projections are described in the text. 



state. We use the ao(980) line-shape parameters from Table IVII Similar details of the 
/o(980) parameterization are unimportant as it is used only in systematic studies. 

For low-mass 7i^7i^{a) and Ktt{k,) S-wave contributions we choose a parameterization 
with a complex pole ttir [T3 |. 

An{m) = (9) 

with m„ = (470 — ^220) MeV and = (710 — z310) MeV, which is adequate for our small 
sample. 

Figure [To] shows the Dalitz plot and three projections for Xci n~^'ir~rj. There are clear 
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contributions from ao(980)^7r^ and /2(1270)?7 intermediate states, and significant accumu- 
lation of events at low vr'^vr" mass. Note that the ao(980) can contribute in two decay 
modes to the Dalitz plot. An isospin Clebsch-Gordan decomposition for this decay, given 
in Appendix |Vll Equation [211 shows that amplitudes and strong phases of both charge- 
conjugated states should be equal. The overall amplitude normalization is arbitrary and we 
set aaQ(98o)+ = aao(980)- = 1- All other fit components are defined relative to this choice. 

Our initial fit to this mode includes only ao(980)^7r''' and /2(1270)r7 contributions, but 
has a vanishing probability of 0.1% for describing the data due to the accumulation of 
events at low tt+tt^ mass. To account for this we try K^, /o(980), p(770), and a resonances. 
Only the p(770) and the a give high fit probability. However the decay xa —>■ p(770)?7 is 
C7-forbidden, and the low- mass distribution is not well represented by the p(770), which 
only gives an acceptable fit due to its large width and the limited statistics of our sample. 
The a describes well the low tt+tt" mass spectrum, and we describe the Dalitz plot with 
ao(980)^7r^, /2(1270)r7, and ar] contributions. Table rvTTl gives the results of this fit, which 
has a probability to match the data of 66%. The angular distributions for ao(980)^ and a 
meson decays are uniform, and for /2(1270) are taken from Table IVl for quantum numbers 
JjL = 121. We assume that a possible contribution from L = 3 is small and it is neglected. 
The systematic uncertainties shown in the table were obtained from variations to this 

TABLE VI: Resonance parameters in the Xci ??7r+7r^ mode comparing their nominal value to 
fit values when individual resonance parameters are allowed to float. 



Parameter 


Nominal Value 


When Floating 


m(ao(980)), MeV/c^ 


999 


1002±18 


grj^, MeV/c2 


620 


637±49 


g^K, MeV/c2 


500 


523±154 


Re(m^), MeV/c^ 


470 


511±28 


Im(m^), MeV/c^ 


-220 


-102±50 



nominal fit as discussed below. We allow the 2D-efficiency to vary with its polynomial 
coefficients constrained by the results of the fit to the simulated events; the mass of the 
ao(980) and its coupling constants are allowed to fioat, the parameters of the cr-pole are 
allowed to fioat, and we allow additional contributions from p{770)ri, /o(980)?7, K^r], and 
vTi (1400)7?. The results of allowing the resonance parameters to fioat as compared to their 
fixed values used in the nominal fit are shown in Table I VII For the additional contributions 
we do not observe amplitudes that are significant and we limit their individual fit fractions 

TABLE VII: Fit results for Xci — * r/TT+vr" Dalitz plot analysis. The uncertainties are statistical 
and systematic. Allowing for interference among the resonances changes the fit fractions by as 
much as 20% in absolute terms as discussed in the text. 



Mode 




Fit Fraction (%) 


ao(980)=^7r=F 


1 


75.1 ±3.5 ±4.3 


/2(1270)7? 


0.103 ±0.014 ±0.005 


14.4 ±3.1 ± 1.9 


ar] 


0.41 ±0.05 ±0.10 


10.5 ± 2.4 ± 1.2 
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FIG. 11: (a) Dalitz plot, (b)-(d) projections on the three mass squared combinations for Xci 
K^K^TT^. The displayed fit projections are described in the text. The contribution from 7r^K~ 
resonances are not shown. They look similar to the ir^K^ resonance components for the relevant 
projections. 



to roughly less than 5% of the Dalitz plot. We note that with higher statistics this mode 
may offer one of the best measurements of the parameters of the ao(980). These results are 
consistent with the substructure analysis of this mode by BES 

The Dalitz plot for Xci K^K~tt^ decay and its projections are shown in Figure [11], 
and for Xci tx^K~K^ in Figure [121 We do a combined Dalitz plot analysis to these modes 
taking advantage of isospin symmetry. An isospin Clebsch-Gordan decomposition for these 
decays, described in Appendix rvT| shows that these two Dalitz plots should have the same 
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FIG. 12: (a) Dalitz plot, (b)-(d) projections on the three mass squared combinations for Xci 
7r~^K~Kg. The displayed fit projections are described in the text. The tt^K~ resonances are not 
shown. They look similar to the n^Kg resonance components for the relevant projections. 

set of amplitudes for all K*K and ao(980)7r intermediate states. The relative factor —V2 
between the two Dalitz plot amplitudes does not matter due to the individual normalization 
of their probability density functions. In the combined fit to these two Dalitz plots, we 
use the following constraints on the amplitudes: ax* = clk*+ = clr*- = clr*^^ = (^j^'O, and 
0^(980) = cta(980)+ = c^a(980)- = (^a{980)° ■ The Overall amplitude normalization is arbitrary and 
we set ax* (892) = 1- 

Visual inspection shows apparent contributions from i^'*(892)^A''^, K* [892)^ Kg, 
K*{U30)^K^, K*{im)^K°s, ao(980)°7r°, and ao(980)±7r^. It is not clear if the K*{U30) 
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are or and many other Kit and KK resonances can possibly contribute. The angular 
distribution for scalar resonance decays is always taken as uniform. The shape for ii'*(892)i^ 
is taken from Table IVl for quantum numbers JjL = 110. We also modeled this contribution 
with the JjL = 112 angular distribution with the difference taken as a systematic uncer- 
tainty. The shape for i^2(1430)if is taken for quantum numbers JjL = 121, with a possible 
contribution from L = 3 ignored. Our best fit result is shown in Table IVIIIt displaying 
statistical and systematic errors. This fit has a good probability of matching the data, 

TABLE VIII: Results of the combined fits to the Xd ^ K+K~tt'^ and Xd ^ ttKK'^ Dahtz plots. 
Allowing for interference among the resonances changes the fit fractions by as much as 15% in 
absolute terms as discussed in the text. 



Mode 


aR 


Fit Fraction (%) 


K*{892)K 


1 


31.4 ± 2.2 ± 1.7 




3.8 ±0.4 ±0.2 


30.4 ±3.5 ±3.7 


K^{U30)K 


0.44 ± 0.06 ± 0.04 


23.1 ±3.4 ±7.1 


ao(980)7r 


6.1 ±0.6 ±0.6 


15.1 ±2. 7± 1.5 



73%, and agrees with fits done to the separate Dalitz plots not taking advantage of isospin 
symmetry. The addition of a nK contribution does not improve the fit probability and the 
amplitude is statistically significant only at the four standard deviation level. Similar 
behavior is noted for a non-resonant contribution. Systematic uncertainties are evaluated by 
observing deviations from the nominal fit as described above for the Xci n^ii^ri analysis. 
When systematics are taken into account neither the kK nor non-resonant contribution has 
a fit fraction significant at the three standard deviation level. Their individual contributions 
are limited to less than roughly 10% of the Dalitz plot. These observations are consistent 
with those by BES fl in the n+K-R^ mode. 

While this model of the contributions to the Dalitz plots gives a good description of the 
data it is clearly incomplete. There should be interference among the resonances, which 
the simple approach of Equation O does not take into account. To quantify the effect of 
interference we have repeated the Dalitz analyses using a "quasi-coherent" (note that is 
always positively-defined) sum of amplitudes with floating complex weighting factors instead 
of aR. This causes the fit fractions to change by as much as ~20% absolutely in the r/vrvr 
analysis and ~15% absolutely in the KKtt analysis. With more data a study of the Xci 
decay substructure taking into account the effect of Xci polarization, using the orientation 
of the Xci decay plane with respect to the Xci flight direction, in a partial wave analysis 
would be a complete description of these decays. The matrix element amplitudes should 
include the partial wave dependent angular distributions and interference effects among the 
resonances. See Ref. [l5| for appropriate prescriptions. 

IV. SUMMARY 

We have searched for and studied selected three-body hadronic decays of the Xco, Xci, and 
Xc2 produced in radiative decays of the ip{2S) in e+e" collisions observed with the CLEO 
detector. Many of the channels covered in this analyses are observed or limited for the 
first time. Our observations and branching fraction limits are summarized in Table IIVI In 
Xci ~^ TT^n^r] we have studied the resonant substructure using a Dalitz plot analysis simply 
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modeling the resonance contributions as non-interfering amplitudes. Our results are summa- 
rized in Table IVTTl We observe clear signals for the 00(980)^7?^ and /2(1270)?7 intermediate 
states, and a low- mass n~^7r~ enhancement. Similarly for Xci KKtt our results are sum- 
marized in Table IVIIIt assuming, based on isospin, that the TiKKg and K~^K~Tr^ plots are 
identical. We observe K*{892)K, K;{U30)K and likely ao(980)7r contributions in the KKtt 
Dalitz plots. Other conclusions about S-wave contributions are likely model-dependent. 
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VI. APPENDIX: CLEBSCH-GORDAN DECOMPOSITION FOR xa 

In order to constrain amplitudes and phases in Xci decays we use a Clebsch-Gordan 
decomposition of the Xci (the state with |/ = 0, = 0)) for possible isospin subsystems: 

Xci ^ {Ktt)j=,/2K, (10) 

Xci (Kn)j=y2K, (11) 

Xci {KK)j=,n. (12) 

Below we use the Clebsch-Gordan decomposition rules, \J,M) = J2f Cf\mi,m2) f from 
Ref. 0. 

A. Clebsch-Gordan decomposition for K* —>■ Kir decays 

We assume that K* mesons with 1=1/2 form two isodoublets: (A'*'*', K*^) and {K*^ , 
K*~) with {Iz = \ilz = ~\) respectively. The Clebsch-Gordan decomposition rules for 
1 X I isospin states are: 

K*^ = y|vr+ir° - ^jl^r^K^ (13) 
K*' = ^n'K' - f^n-K\ (14) 
r° = ^n^K- - ^vr°Z°, (15) 
K*- = J^rr'K- - f^^'^ ■ (16) 
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B. Cases of Xci (ifvr)j=i/2K and Xd {Ktt)i=i/2K decays 
For Xci K*K and xa K*K modes we use the | x | rule: 



1 

V2' 



x,, = ^{k*'k'-k*-k' 



Combining Equations [T7] and [18] with Equations [T3lfT6l we get 



XC1V2-- 
XC1V2 



n^K-)K° + (7i-K")K 



Assuming charge symmetry the amphtudes in Equations [19] and [20] should be equal, 
these equations we get the ratio of rates: 

Tixci ^ n+K-K^)/T{xci ^ K+R-n'^) = 2, 



-7r°) = 2, 



or their sum 



C. Case of Xci {KK)j=itt decay 



For Xci c-vr modes we use the 1x1 rule: 

1 



Xci 



For a — > KK we use the ^^^^s: 



a+ = K+K 



a' = K^K-. 

Combining Equation [24] with Equations [25l[27] we get 

V2 L 



From Equation [28] we get the ratio of rates: 

Tixci n+K-K^)/T{xci ^ K+R-n') = 2, 



or their sum 



nx 



^K-K") + r(xci ^ TT^ir+iT) = 4 ■ r(xci ^ K+K-n 
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D. Consequences for Dalitz plot analysis 



Comparing Equations [2T1I231 for intermediate states with K* and Equations [291I3T] for 
intermediate states with a(980), we note that they are identical Thus observations of 

n^K'K^ and charge conjugated t[~ K^K^ final states on the same Dahtz plot will yield the 
certain ratio between K* and a(980) amplitudes for the irKK^ Dalitz plot. The same ratio 
between K* and a(980) amplitudes is expected for the K~^K~n^ Dalitz plot. This isospin 
analysis implies that these two Dalitz plots, nKKg and K~^K~n^, can be parametrized using 
a common set of parameters for each K* and a(980) intermediate state. From Equations! 
we can write equations between decay amplitudes with K* mesons 



for K^K 7r° : = ——j^R*+R-^R+R-T,':i = —-^ax'-R+^R+R-nO, (32) 
where the signs assume equal phases 

(33) 

Similar equations between amplitudes with ao(980) can be obtained from Equation [28] 

for TlKK^g : aa(980)+7r-^7r/^E'0 = <^a{980)-T7+^nRR0 = 



)R*+ = (PR*- = (PR*o = 9k*o. 



for K^K 7r° : = ^aa{9SQ)0n0^R+R-7r0, (34) 



(35) 



1 

assuming equal phases 

0a{98O)+ = 0a(98O)- = 0a(98O)O- 

Equations [32] and [34] predict that the ratio of amplitudes between these two Dalitz plots is 
— V^. This relative factor does not matter, because each Dalitz plot is normalized separately. 
In the combined fit we ignore this — \/2 factor between the amplitudes in the nKKg and 
K^K'TT^ Dalitz plots, and use common fit parameters for each of K* mesons and ao(980) 
intermediate states 

ax* = clr*+ = clk*- — o,R*o = CL-^to, (36) 

Cta(980) = 0-0(980)+ = '^0(980)- = '2a(980)0 • (37) 

In Sections IVI CI and IVIBI we have checked an isospin symmetry between ttKK^ and 
K^K~TT^ decays for two particular intermediate states. However, this symmetry is valid 
independently of how the amplitudes are decomposed into two-body sub-amplitudes, and 
thus should be valid for all points on the Dalitz plot. 
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